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ABSTRACT 

The  mixed  mode,  near-field  state  of  stresses  surrounding  a  crack 
propagating  at  constant  velocity  is  used  to  derive  a  relation  between  the 
dynamic  stress  intensity  factors  Kj,  Kjj,  the  remote  stress  component  oqx  and 
the  dynamic  isochromatics.  This  relation  together  with  an  overdetermini stic 
least-square  method  form  the  basis  of  a  data  reduction  procedure  for  extract¬ 
ing  dynamic  KI*  KII  and  aQx  from  the  recorded  dynamic  photoelastic  pattern 
surrounding  a  running  crack.  The  overdeterministic  least-square  method  is 
also  used  to  fit  static  isochromatics  to  the  numerically  generated  dynamic 
isochromatics.  The  resultant  static  Kj,  Kjp  and  aQX  are  compared  with  the 
corresponding  dynamic  values  and  estimates  of  errors  involved  in  using  static 
analysis  to  process  dynamic  isochromatic  data  are  obtained.  The  data  reduc¬ 
tion  procedure  is  then  used  to  evaluate  the  branching  stress  intensity 
factor  associated  with  crack  branching  and  the  mixed  mode  stress  intensity 
factors  associated  with  crack  curving. 

INTRODUCTION 

In  recent  years,  two  dimensional  dynamic  photoelasticity  has  been  used 
to  determine  experimentally  the  dynamic  stress  intensity  factor  surrounding 
a  propagating  crack  and  to  establish  a  dynamic  fracture  toughness,  Kjq,  ver¬ 
sus  crack  velocity,  a,  relation  which  is  believed  to  control  dynamic  fracture. 
An  excellent  article  on  this  use  of  dynamic  photoelasticity  for  studying 
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dynamic  fracture  has  been  written  recently  by  J.  W.  Dally  [1]*.  Since  fracture 
resistance,  i.e.  dynamic  fracture  toughness  «ID  is  equal  to  the  driving  force, 
i.e.  instantaneous  dynamic  stress  intensity  factor,  during  crack  propagation 
the  dynamic  stress  intensity  facto-',  which  is  extracted  from  transient  dynamic 
isochromatics  surrounding  the  propagating  crack  tip,  is  used  to  measure  the 
dynamic  fracture  toughness.  A  commonly  used  data  reduction  procedure  for 
this  purpose  is  to  fit  a  theoretical,  near-field,  static  isochromatics  to  the 
recorded  experimental  dynamic  isochromatics  and  to  then  equate  the  resultant 
static  stress  intensity  factor  of  the  former  to  the  unknown  dynamic  stress 
intensity  factor  of  the  latter  [2-5].  Error  estimates  for  using  a  static 
near-field  stress  to  extract  the  dynamic  stress  intensity  factor  have  been 
made  by  several  investigators  [6-8]  and  in  particular,  exhaustively  by 
Rossmanith  and  Irwin  [8]. 

Studies  of  the  static  isochromatic  patterns  under  mixed  mode  loading  con¬ 
ditions,  i.e.  in  the  presence  of  combined  Kj.  and  crack  tip  deformation, 
were  made  by  Smith  and  Smith  [9],  Gdousto  and  Theocaris  [10]  and  more  recently 
by  Dally  and  Sanford  [11,12],  These  isochromatics  are  all  characterized  by 
their  unsymmetric  patterns  with  respect  to  the  crack  line.  While  the  simple 
static  isochromatics  used  to  demonstrate  the  importance  of  mixed  node  loading 
conditions  in  the  above  references  were  mathematically  Contrived,  Klein  deter¬ 
mined  Kj  and  Kjj  from  actual  photoelastic  patterns  of  a  crack  approaching  a 
hole  [13].  Literature,  however,  is  scant  on  other  mixed  mode  photoelastic 
investigations  despite  the  well  studied  theoretical  conditions  under  which 
such  unsymmetric  isochromatics  can  exist. 

On  the  other  hand,  dynamic  isochromatics  surrounding  a  running  crack 

often  exhibits  moderate  unsymmetry  but  such  photoelastic  patterns  were 
* 

Number  in  brackets  refer  to  Reference  at  the  end  of  this  paper. 


heretofore  considered  experimental  abnormalities  and  were  ignored  by  averaging 
the  unsymmetric  to  symmetric  patterns  during  the  data  reduction  process.  Care¬ 
ful  postmortem  inspection  of  the  fracture  specimens,  however,  show  that  the 
slightly  unsymmetric  isochromatics  are  often  associated  with  slightly  curved 
crack  propagation  paths  which  undoubtedly  are  caused  by  the  small  dynamic  Kjj, 
coexisting  with  the  dominating  dynamic  Kj  value.  This  effects  is  akin  to  the 
small  but  noticeable  influence  of  a  small  Kjj  on  fatigue  crack  propagation 
reported  fifteen  years  ago  [14].  Grossly  unsymmetric  isochromatics  associated 
with  large  crack  curving,  crack  branching  and  propagating  multiple  cracks,  on 
the  other  hand,  are  unmistakably  caused  by  larger  values.  The  exact  rela¬ 
tion  between  the  amount  of  crack  curving  and  the  dynamic  Kj  and  associated 
with  the  propagating  crack  tip  would  provide  a  dynamic  crack  propagation 
law  under  mixed-mode  crack  tip  deformation  similar  to  the  KjD  versus  8  rela¬ 
tion  under  consideration  for  Mode  I  dynamic  crack  propagation  [15],  Since 
the  relative  magnitude  of  dynamic  usually  is  small  with  respect  to  the 
dynamic  Kj  values  for  many  crack  propagation  problems,  these  values  must 
be  determined  accurately  if  a  meaningful  dynamic  mixed-mode  crack  propagation 
law  is  to  be  deduced  from  these  results.  The  purpose  of  this  investigation 
is  to  develop  a  data  reduction  procedure  for  extracting  Kj  and  from 
recorded  dynamic  isochromatics  surrounding  a  running  crack  tip.  The  developed 
procedure  is  also  used  to  estimate  errors  involved  in  using  static  mixed-mode 
crack  tip  stress  field  in  place  of  the  corresponding  dynamic  field. 

THEORY 

The  three  stress  components,  o  ,  o  ,  and  0  in  terms  of  the  Modes  I 

xx  yy  xy 

and  II  dynamic  stress  intensity  factors,  Kj  and  K^,  plus  the  remote  stress 
★ 

component,  o  ,  the  maximum  shear  stress  can  be  represented  as  [16,17,18] 

U  A 

★ 

Note  that  the  sign  conventions  of  o„v  is  negative  of  that  used  in  References 
[2,7,8,11,12].  ox 
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r  and  e  are  the  traditional  polar-coordinates  with  origin  at  the  moving  crack 
tip  and  c,  c^ ,  and  c2  are  the  crack  velocity,  dilatational  and  distortional 
wave  velocities,  respectively. 

For  each  given  e  the  radial  distance,  r,  can  be  explicitly  solved  from 
equation  (la)  for  a  given  crack  velocity,  c,  mixed  mode  dynamic  stress  inten¬ 
sity  factors  of  Kj  and  K^,  and  remote  stress  component  of  oQX.  Thus  the 
theoretical  isochromatics  xm  for  a  straight  crack  propagating  under  constant 
velocity  can  be  easily  constructed.  When  equation  (1)  is  used  to  extract  Kj 
and  Kjj  from  recorded  dynamic  isochromatics,  one  can  generate  9000  plus 
isochromatic  loops  for  various  combinations  of  Kj,  K^,  and  oQX  following 
Etheridge  et  al .  [19]  and  use  a  search  routine  to  match  the  computer  stored 
isochromatics  with  the  experimental  isochromatics.  In  this  paper,  however,  the 
direct  and  overdeterministic  procedure,  which  was  developed  by  Bradley  [3] 
and  updated  with  improved  numerical  techniques  by  Sanford  et  al .  [12],  of  least 
square  fitting  equation  (1)  to  a  large  number  of  measured  isochromatic  data 
point  is  used. 

For  the  above  overdeterministic  method  of  least  square  fitting  a  theoreti¬ 
cal  isochromatics  im  at  coordinate  locations  (r^,  e^),  equation  (1)  can  be 
rewritten  in  a  functional  form,  F^,  which  represents  the  deviation  between 
the  square  of  the  calculated  and  measured  isochromatics 
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where  k  =  1,2 . .  is  the  number  of  data  points  on  the  recorded 

isochromatic. 
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Obviously  =  0  when  the  theoretical  and  measured  isochromatic  coincides. 
In  general,  such  coincidence  cannot  be  expected  and  thus  Kj,  Kj j  and  oQX 
must  be  optimized  in  order  to  minimize  the  total  summation  of  F^.  For  such 
minimization,  consider: 
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The  subsequent  steps  of  extracting  the  optimized  Kj ,  Kjj,  and  aQX  using 
the  least  square  method  is  a  standard  numerical  procedure  which  is  also  iden¬ 
tical  to  that  of  Reference  [12]  and  thus  will  not  be  repeated.  Basically, 
the  procedure  consists  of  a  sequence  of  error  reduction  estimates  aKj,  aKjj, 
and  Aa  starting  from  an  initially  estimated  KT,  KTT,  and  o  .  The  numerical 
procedure  has  considered  to  converge  when  the  current  AKj,  AKj j ,  and  Aoqx  are 

less  than  0.1%  of  the  previously  calculated  KT,  KTT  and  o  . 

Ill  ox 

THEORETICAL  ISOCHROMATICS 

As  one  checks  on  the  accuracy  of  equation  (1),  the  dynamic  isochromatics 
should  coincide  with  that  of  Reference  [11]  in  the  limiting  case  of  crack 
velocity  c+0.  Also  informative  would  be  a  comparison  between  this  static 
isochromatics  and  the  corresponding  dynamic  isochromatics  for  the  identical 
Kj,  Kjj,  and  oqx  for  a  given  crack  velocity. 

In  order  to  provide  such  theoretical  comparison,  theoretical  isochromatics 
were  generated  using  the  same  model  fringe  constant  of  fQ/h  =  1.73  MPa/fringe 
(250  psi/fringe),  Kj  =  .879  MPa  y/m  (800  psi  /in’)  and  Kjj  of  Reference  [12]. 
Static  isochromatics  were  approximated  by  setting  c  =  0.0001  Cj  in  equation 
(1).  A  constant  crack  velocity  of  c  =  0.15  Cj  was  chosen  to  generate  the 
dynamic  isochromatics.  In  addition  Cj  =  2,400  m/sec  (94,300  in/sec)  and 
c2  =  1,160  m/sec  (45,800  in/sec)  were  assumed  in  the  dynamic  analysis.  In 
the  following  some  representative  static  and  dynamic  isochromatics,  which 
were  originally  plotted  on  a  CALCOM  plotter  and  then  inked  for  presentation, 
are  shown. 

Mode  I  static  and  dynamic  isochromatics  have  been  studied  by  others 
[6,8,20]  and  thus  many  of  the  sample  calculations  in  these  references  were 
reproduced  as  a  partial  check  of  the  computer  code  developed  for  plotting 
dynamic  isochromatics.  In  particular,  the  dynamic  isochromatics  generated 
in  Reference  [20]  were  all  reproduced  with  the  new  algorithm  and  the  two 
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results  agreed  within  0.01  percent.  For  comparison  purpose,  some  of  the 
results  of  pure  Mode  I,  pure  Mode  II  and  mixed-mode  static  isochromatics  in 
Reference  [12]  together  with  their  dynamic  counterpart  are  listed  in  the 
following. 

Figures  1,  2  and  3  show  typical  static  and  dynamic  isochromatics  of  fringe 
orders,  2.5  and  3.5  surrounding  a  crack  tip.  As  amply  discussed  by  others 
[7,  8,  11],  for  negative  values  of  oqx  the  isochromatics  lean  forward  in  the 
direction  of  crack  propagation  while  for  positive  values  of  oQx  the  isochrom¬ 
atics  lean  backward  away  from  the  direction  of  crack  propagation.  The  static 
isochromatics  in  Figure  1  coincides  with  that  in  Figure  4  of  Reference  [11]. 
Figures  2  and  3  show  further  backwards  leaning  isochromatics  which  are  due  to 
the  higher  oqx/Kj  =  0.5  and  0.75.  As  shown  in  this  figure,  the  differences 
between  the  static  and  dynamic  isochromatics  are  small  for  forward  leaning 
loops  in  Figure  1  but  increase  slightly  with  increase  in  backward  leaning 
angles  in  Figures  2  and  3. 

Figure  4  shows  the  static  and  dynamic  isochromatics  of  N  =  2.5  and  3.5 
for  a  pure  Mode  II  crack  tip  deformation.  The  static  results  are  in  agree¬ 
ment  with  that  in  Figure  5  of  Reference  [11].  Unlike  the  pure  Mode  I  crack 
tip  deformation,  little  difference  between  the  static  and  dynamic  isochromatics 
is  observed.  Experimentally,  a  crack  has  seldom  been  observed  to  propagate 

in  a  straight  line  under  pure  Mode  II  crack  tip  deformation  but  will  generally 

* 

curve  to  achieve  a  maximum  Kj  with  attendant  smaller  K^.  Thus  a  pure  Mode  II 
dynamic  isochromatic  may  be  a  mute  academic  study.  Pure  Mode  II  isochromatics 
of  a  stationary  crack  were  reported  in  Reference  [21]. 


Since  the  near-field  dynamic  stress  field  used  to  generate  the  dynamic 
isochromatics  of  equation  (la)  is  for  a  straight  crack,  the  extent  of  near¬ 
field  for  a  curved  crack  must  be  reduced  accordingly  for  this  straight  crack 
approximation  to  hold. 
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Probably  the  most  frequently  encountered  mixed-mode  static  and  dynamic 
isochromatics  are  generated  by  the  presence  of  relatively  small  Kjj  and  0qx 
coexisting  with  the  Kj  mode  of  crack  tip  deformation.  Figure  5  shows  a  typical 
asymmetric  isochromatic  patterns  where  some  differences  between  the  static  and 
dynamics  isochromatics  are  noted.  Figure  6  shows  significant  change  in  iso¬ 
chromatic  patterns  when  the  sign  of  the  remote  stress  component  oQx  is  changed. 
Needless  to  say  the  static  isochromatics  coincide  with  those  in  Figures  8a  and 
8c,  respectively  in  Reference  [11]. 

Other  static  isochromatic  patterns  shown  in  Reference  [11]  together  with 
the  corresponding  dynamic  isochromatics  were  also  generated  mainly  to  verify 
the  computer  plotted  isochromatics  for  c  =  0.0001  c-j  cases.  Details  of  the 
changes  in  the  shapes  of  the  static  isochromatics  are  discussed  thoroughly  in 
Reference  [11]  and  thus  will  not  be  reproduced  here.  In  essence  the  correspond¬ 
ing  dynamic  isochromatics  are  found  to  be  slightly  larger  and  essentially  follow 
the  general  shape  of  the  static  isochromatics  similar  to  results  shown  in 
Figures  1  through  6. 

ERROR  ANALYSIS 

As  mentioned  previously,  the  basic  data  reduction  procedure  developed 

in  this  investigation  consists  of  least-square-fitting  the  mixed  mode, 

dynamic  near-field  state  of  stress  to  a  large  number  of  data  points  of  a 

dynamic  isochromatic  field.  For  a  known  constant  crack  velocity  then,  this 

overdetermini stic  numerical  procedure  provided  the  associated  dynamic  Kj, 

Kjp  and  aQx.  Internal  consistency  of  this  procedure  was  first  verified  by 

least-square  fitting  a  dynamic  isochromatic  every  =  18,36 . ,360,  to 

the  theoretical  dynamic  isochromatic  generated  for  given  Kj ,  Kjj,  and  aQx 

and  then  recovering  these  three  dynamic  parameters  through  the  overdetermin- 

istic  method.  The  recovered  KT,  KTT  and  a  werewithin  0A%  of  their 

ill  ox 

original  values  for  all  five  cases  represented  by  Figures  1  through  6. 
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The  differences  in  static  and  dynamic  isochromatics  for  the  commonly 
observed  maximum  crack  velocity  of  about  c  =  0.15  c^ ,  as  shown  in  Figures  1 
through  6,  are  relatively  small.  Nevertheless,  a  quantitative  assessment  of 
the  error  involved  in  using  static  analysis  to  extract  dynamic  Kj  and  is 
necessary  since  static  analysis  has  been  used  extensively  for  reducing  the 
dynamic  photoelasticity  data  associated  with  a  running  crack. 

Errors  involved  in  using  a  three-parameter  static  near-field  state  of 
stress  to  characterize  a  dynamic  isochromatic  field  associated  with  a  constant 
velocity  crack  were  estimated  by  least  square  fitting  a  stt.;c  isochromatics 
(i.e.  c  =  0.0001  c-|)  to  theoretical  dynamic  isochromatics  wh.ch  were  generated 
through  the  use  of  equation  (1).  The  convnonly  used  data  reduction  procedure 
was  simulated  by  least  square  fitting  the  static  field  at;  1)  three  points 
of  em„„*and  0m,„  ±  6°  along  the  two  dynamic  isochromatic  lobes  on  both  sides 
of  the  crack  and  2)  along  10  points  each  straddling  0  and  between  the 
maximum  width  of  the  two  isochromatic  lobes  on  both  sides  of  the  arch.  The 
difference  between  the  resultant  static  Kj,  Kjj,  and  oQX  and  the  corresponding 
dynamic  values  then  constitute  the  theoretical  errors  involved  in  using  the 
static  near-field  isochromatics  for  data  reduction. 

Figure  7  shows  the  progressive  increase  in  error  with  increased  crack 
velocity  for  pure  Mode  I  crack  tip  deformation.  Only  the  20  point,  i.e.  10 
points  each  on  each  of  the  two  symmetric  isochromatic  lobes,  data  fitting 
procedure  was  used  for  obtaining  these  results  as  well  as  those  in  Figure  7. 
Figure  7  shows  that  errors  involved  in  the  overdeterministic  method  of  least- 
square  fitting  the  near-field  static  isochromatics  of  two  parameters,  Kj  and 
°ox  with  KII  =  °*  does  not  differ  substantially  with  oQX,  with  the  exception 

★ 

For  definition  of  e  , 

FTlaX 


see  References  [2, 3, 5, 7, 8]. 
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of  o0X/Kj  =  ±  0.75,  in  the  crack  velocity  range  of  c/Cj  <  0.2  and  this  finding 
is  in  agreement  with  that  of  Reference  [7].  The  lack  of  systematic  changes  in 
errors ,  i  .e.  for  aQX/Kj  =  ±  0.75  is  an  indication  that  the  optimization  process 
involved  in  the  least-square  method  is  functioning  properly.  Significant,  how¬ 
ever,  is  the  systematic  incresase  in  the  error  amounting  to  13-24  percent  of 
the  correct  dynamic  values  at  extremely  high  crack  velocity  of  c  =  0.25  c^ . 

For  a  crack  velocity  of  c/Cj  =  0.15  the  estimated  errors  shown  in  Figure  7 
are  within  2^4  percent  of  those  of  Rossmanith  and  Irwin  [8]  who  conducted  a 
more  extensive  error  analysis  involved  with  Mode  I  dynamic  crack  propagation. 

In  the  unlikely  event  of  a  dynamic  crack  propagation  under  pure  Mode  II 
crack  tip  deformation,  a  two  parameter  static  stress  field,  Kjj  and  oQx  with 
Kj  =  0,  can  be  fitted  to  this  dynamic  isochromatics.  The  errors  involved  in 
the  fitting  process  is  shown  in  Figure  8.  Unlike  the  pure  Mode  I  crack  tip 
deformation,  the  errors  involved  in  using  the  static  near-field  stress  are 
small  as  may  have  surmised  from  Figure  4. 

Errors  estimation  in  fitting  a  static  near-field  mixed  mode  isochromatic 
to  a  mixed-mode  dynamic  isochromatic  will  obviously  vary  with  the  relative 
magnitudes  of  Kj ,  Kjj  and  oqx.  In  order  to  estimate  such  errors,  the  static 
three  parameter  stress  field  was  least-square  fitted  to  an  arbitrary  mixed¬ 
mode,  dynamic  isochromatic  associated  with  Kj  =  .879  MPa  ^  (800  psi  /Tn) 
and  Kjj  =  .219  MPa  </m  (200  psi  /Tn)  with  varying  oQX.  Figures  9a  and  9b  show 
the  widely  varying  errors  in  the  Kj  and  Kjj  determined  by  such  optimization 
procedure.  Of  particular  concern  is  the  unpredicted  large  errors  in  Kj  and 
associated  with  higher  crack  velocities.  While  this  large  error  could 
in  part  be  due  to  the  distorted  shape  of  isochromatics  as  shown  in  Figure  6, 


* 


These  values  were  selected  from  Reference  [11]. 


these  results  nevertheless  indicate  the  importance  of  using  the  dynamic  near 
field  stress  for  extracting  dynamic  Kj  and  Kjj  when  the  isochromatics  are 
distorted. 

EXAMPLES 

The  developed  computer  code  for  computing  dynamic  Kj,  K^,  and  oQx  was 
then  used  to  determine  the  dynamic  Kj  and  from  the  dynamic  isochromatics 
associated  with  a  curving  crack  and  a  branching  crack. 

Figure  10  shows  two  frames  out  of  a  16-frame  dynamic  photoelastic  record 
of  a  curving  crack  impacting  a  12.7  mm  (1/2-inch)  diameter  hole  in  a  9.53  mm 
(3/8-inch)  thick,  254  x  254  mm  (10  x  10  inch)  Homalite-100  plate  loaded  under 
fixed  gripped  tension.  The  crack  emanated  from  a  small  precrack  at  the  top 
edge  of  plate  upon  reaching  a  critical  load  and  propagated  slightly  off  the 
centerline  until  it  was  pulled  into  the  higher  tension  field  surrounding  the 
hole.  Further  details  of  the  experimental  stepup, crack  velocity  measurements 
and  dynamic  calibration  of  the  Homalite-100  material  used  are  found  in 
Reference  [22]. 

Figure  11  shows  three  frames  out  of  a  16-frame  dynamic  photoelastic  record 
of  a  crack  propagating  and  branching  in  a  3.18  mm  (1/8-inch)  thick,  254  x  254  mm 
(10  x  10  inch)  Homalite-100  plate  loaded  under  fixed-gripped  tension.  Again 
details  of  the  experiment  can  be  found  in  Reference  [22]. 

Figure  12  shows  the  dynamic  Kj  and  Kjj  variations  obtained  from  the 
four  dynamic  photoelastic  patterns  preceding  crack  curving  and  crack  impacting 
the  hole  as  shown  in  Figure  10.  Although  not  conclusive,  the  ratio  of  dynamic 
Kjj/Kj  increases  rapidly  prior  to  crack  curving,  thus  leading  to  a  natural 
speculation  that  a  moderately  mixed  mode  local  dynamic  state  of  stress 
results  in  curving  of  a  propagating  crack. 
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Figure  13  shows  the  dynamic  Kj  and  Kjj  for  three  branches  of  the  cracks 
shown  in  Figure  11.  While  the  continuous  "right  branch"  crack  shows  moderate 
changes  in  dynamic  Kj  and  K^,  the  left  branch  shows  a  sharp  drop  in  dynamic 
Kj  after  crack  branching.  By  extrapolating  the  dynamic  KjS  associated  with 
"left  branches  Nos.  1  and  2",  an  after-branching  dynamic  stress  intensity 
factor  of  0.923  MPa  v'riT  (840  psi  /Tn)  is  obtained.  The  branching  stress  inten¬ 
sity  factor,  i.e.  immediately  prior  to  branching,  is  estimated  to  be 
2.03  MPa  ►¥  (1850  psi  /in).  Although  this  singular  data  is  higher  than  the 
branching  stress  intensity  factor  of  1.38  MPa  /m  (1250  psi  /Tn)  quoted  in  Refer¬ 
ence  [1],  the  ratio  of  2.03/0.923  =  2.2  is  consistent  with  the  postulate  crack 
branching  occurs  to  dissipate  fracture  energy  along  two  propagating  cracks. 

It  is  also  interesting  to  note  that  dynamic  which  is  a  relatively  low 
0.11  MPa  Jm  (100  psi  /Tn)  prior  to  crack  branching  nearly  doubles  immediately 
after  crack  branching  and  is  consistent  with  the  static  results  of  Reference 
[23]. 

CONCLUSIONS 

1.  A  data  reduction  scheme  for  evaluating  the  unsymmetric  dynamic 
isochromatic  associated  with  dynamic  mixed-mode  crack  propagation  has  been 
developed. 

2.  Errors  involved  in  using  a  mixed-mode  static  near-field  stress  to 
evaluate  dynamic  Kj  and  from  a  dynamic  unsymmetric  isochromatics  could 
be  larger  than  that  involved  for  pure  Mode  I  or  pure  Mode  II  crack  tip 
deformation. 

3.  The  utility  of  the  developed  data  reduction  procedure  was  demonstrated. 
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FIGURE  I.  STATIC  AND  DYNAMIC  ISOCHROMATICS  NEAR  A 
STATIONARY  (C  =  O.OOOIC,)  AND  PROPAGATING 
( C -  0. 1 5 C | )  CRACK  TIP. 

KjsSOOpsi  yin,  Kns  0,  crox  s  -  0.25  Kj  psi 


FIGURE  2.  STATIC  AND  DYNAMIC  ISOCHROMATICS  NEAR  A 
STATIONARY  (CBO.OOOIC| )  AND  PROPAGATING 
(C*O.I5C,)  CRACK  TIP. 

KjsSOOpsi^A”  t  Kn*0,  <^x  8  O.SKjpsi . 


FIGURE  3.  STATIC  AND  DYNAMIC  ISOCHROMATICS  NEAR  A 
STATIONARY  (C»O.OOOIC,)  AND  PROPAGATING 
IC*0.  I5C,)  CRACK  TIP. 

KT*  800  psi  yin*  Kn»0  ,  s  0.75  K  psi. 


FIGURE  4.  STATIC  AND  DYNAMIC  ISOCHROMATICS  NEAR  A 
STATIONARY  (C  *0.0001  C.)  AND  PROPAGATING 
(C  *0.I5C,)  CRACK  TIP. 

Kj  =  0,  Kns  800psi -/In,  crox* -0.25 Kn  psi 


C  =  O.OOOIC, 
C s0. 15C, 


FIGURE  5.  STATIC  AND  DYNAMIC  ISOCHROMATICS  NEAR  A 
STATIONARY  (C  =  O.OOOI  C.)  AND  PROPAGATING 
(C  =  O.I5C,)  CRACK  TIP. 

Kjs800psi«/»n*  Kns0.25Kipsi-/ln,  o-OJt= -0.25Kj  psi 


FIGURE  6.  STATIC  AND  DYNAMIC  ISOCHROMATICS  NEAR 
A  STATIONARY  (C*  0.000  1C,)  AND  PROPAGATING 
(Cs  0.15  C,)  CRACK  TIP. 

Kj  8  800psi>/firi,  Kn80.25  Kj  psi-v/Tn,  cr0)l8  -0.25  Kj  psi 
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FIGURE  8.  ERROR  IN  DYNAMIC  DETERMINATION  USING  A 
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FIGURE  9 a. ERRORS  IN  DYNAMIC  ^DETERMINATION  USING  A  STATIC 
THREE  PARAMETER  K*  -  K n  -  cfcx  ,  STRESS  FIELD. 
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FIGURE  9b  .  ERRORS  IN  DYNAMIC  DETERMINATION  USING  A 

STATIC  THREE  PARAMETER  KT-Kn-orox  .STRESS  FIELD. 
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FIGURE  10  .TYPICAL  DYNAMIC  PHOTOELASYSf  PATTERNS  IN 
HOMALITF- 100  SINGLE  EDGED  MOT  cm  SPECIMEN 
(FIXED  SHIP  LOADING)  CRACK  A  v  POACHING  A 
CENTRA;,  HOLE,  NO.  W0902’Y). 


(a)  EIGHTH  FRAME  138  /x  SECONDS 


(b)  ELEVEN!  H  FRAME  225  fj.  SECONDS 


FIGURE  II.  TYPICAL  CRACK  BRANCHING  DYNAMIC  PHOTOELASTlC 
PATTERNS  IN  HOMALITEIOO  SINGLE  EDGEDNOTCH 
SPECIMEN  (FIXED  GRIP  LOADING),  NO.  W082270, 
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